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We present a new description of discrete space-time in 1+1 dimensions in terms of a set of 
elementary geometrical units that represent its independent classical degrees of freedom. This is 
achieved by means of a binary encoding that is ergodic in the class of space-time manifolds respecting 
coordinate invariance of general relativity. Space-time fluctuations can be represented in a classical 
lattice gas model whose Boltzmann weights are constructed with the discretized form of the Einstein- 
Hilbert action. Within this framework, it is possible to compute basic quantities such as the Ricci 
curvature tensor and the Einstein equations, and to evaluate the path integral of discrete gravity. 
The description as a lattice gas model also provides a novel way of quantization and, at the same 
time, to quantum simulation of fluctuating space-time. 

PACS numbers: 04.60.-m, 45, 89.70.-a, 31.15.xk 



A central task in defining a discrete version of space- 
time is to identify the relevant degrees of freedom. Addi- 
tionally, any theory of discrete spacetime must take into 
account coordinate invariance, which is a fundamental 
property of General Relativity (GR). This renders the 
identification of fundamental degrees of freedom non- 
trivial. Solving this problem corresponds to fixing the 
gauge of the coordinate invariance symmetry of GR at 
the discrete level. 

Here we address this question and present a binary 
description of space-time in 1+1 dimensions. To do so, 
we borrow the discretization method of causal dynami- 
cal triangulation (CDT) [1]. Following this route leads 
to a formulation of spacetimes in terms of a statistical 
mechanical model, which we identify with a lattice gas 
model [2]. Fluctuations of spacetime are thereby rep- 
resented by different states of the lattice gas, and the 
Boltzmann factor for this statistical model is given by 
the discretized action of general relativity [3] . 

The central idea of this work is to replace the dynam- 
ical geometry of space-time with a fixed physical lattice 
with binary degrees of freedom at its vertices. This con- 
struction allows us to digitalize the geometrical informa- 
tion content of space-times. More precisely, we show how 
to construct a foliated triangulation T from a bit array 
A by using 'forks', and conversely how to construct a 
bit array A from a foliated triangulation T. For the lat- 
ter construction, one first maps T to its dual T*, T* to 
an integer string S, and S to A, as indicated in the fol- 
lowing commuting diagram (the details of which will be 
explained below). 



Triangulation T 



Dual triangulation T* 



forks 



Bit array A 



Integer string S 



(1) 



In addition, we will translate the Pachner moves to op- 
erations on the integer string encoding. (The Pachner 
moves are transformations in the set of triangulations 
which are ergodic, that is, any two triangulations are re- 
lated by a finite sequence of Pachner moves [5].) Finally, 
we will use our binary encoding to formulate meaningful 
quantities of discrete gravity in 1+1 dimensions in the 
natural language of information processing. 

In order to establish these results, let us first recall 
some basic properties of CDT [1]. In CDT the continu- 
ous manifold of space-time is approximated by a piece- 
wise linear manifold [6], where the edge lengths of the 
simplices are assumed to be constant (space-like edges 
have length Z 2 dgo and time-like edges — c^ 2 dgc , we shall 
henceforth assume a = — 1, which is one of the points 
in the 'Euclidean sector' that gives rise to an interesting 
new phase [4]). Additionally, only manifolds that obey a 
global proper discrete time are considered, i.e. manifolds 
with a discrete global time foliation. On the simplicial 
manifold one can define curvature, an action, and other 
quantities [6]. 

We shall here focus on the case of 1+1 dimensions, in 
which the configuration space 7t is formed by all foliated 
triangulations T with with t discrete proper-time steps 
(see Fig. 1). Additionally, we restrict ourselves to sim- 
plicial manifolds with a fixed topology, such as S* 1 x S 1 
(periodic boundary conditions (pbc) in time and space), 
or [0, 1] x S 1 (open boundary conditions in time and pbc 
in space). We consider only connected simplicial man- 
ifolds; this implies, in particular, that there is at least 
one simplex per spatial slice. The topology of the simpli- 
cial manifold is characterized by the Eulcr characteristic 
X ■= N — Ni + N 2 , where N , Ni,N% is the number of 
vertices, edges and faces of T, respectively. Note that in 
2 dimensions x is related to the curvature via the Gauss- 
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FIG. 1: CDT considers only triangulations that obey a global 
time foliation, (a) A foliated triangulation corresponding to a 
flat geometry, since the coordination number of every vertex 
is 6, excluding the boundaries (see Eq. (2)). Any other foli- 
ated triangulation, such as (b), can be obtained by applying 
a sequence of the Pachner moves Rule 1 and Rule 2. 



Bonnet theorem. 

In this setting, the Pachner moves consist of two op- 
erations (called Rule 1 and Rule 2; see Fig. 1). The 
curvature Ri is evaluated at every vertex i: 



R, 



6 



(2) 



where Ci is the coordination number (i.e. number of ad- 
jacent vertices) of i. The discretized Einstein-Hilbert 
action takes the form 



S(T) = 7X - kN 2 , 



(3) 



where 7 and k are coupling constants related to Newton's 
constant Gn and the effective cosmological constant A 
in the continuum, respectively [1]. The quantization is 
carried out by means of a path integral formulation of the 
action (3), which (for a specific topology such as [0, 1] x 
S 1 ) takes the form 



E 



1 



C(T) 



(4) 



Here C (T) is the order of the automorphism group of T, 
and is the remnant of coordinate invariance of GR in the 
discrete theory. 

Considerable progress has been made over the last few 
years — in terms of analytic in I + I [7] and numerical 
studies in 2+1 [8] and 3+1 [4] dimensions — to investigate 
the continuum limit of CDT. Interesting connections with 
other continuum quantum gravity approaches have also 
been established [9-11]. A connection between foliated 
triangulations in terms of random walks was established 
in [12]. Finally, let us mention that other discrete models 
of quantum gravity have been proposed, e.g. [13, 14]. 



From bit arrays to triangulations. We now show how 
a bit array A encodes a foliated triangulation T, as indi- 
cated in Diagram (1). The key observation is that every 
foliated triangulation T is built entirely (except for the 
boundaries) out of certain building blocks that we call 
'forks'. The idea is the following: while, obviously, the 
building block of a general triangulation is the triangle, 
the basic unit of a foliated triangulation is a pair of tri- 
angles that share a space-like edge. We identify the fork 
with this unit; more precisely, each fork consists on 1 
'center', 3 legs and 2 faces (see little diagram on bottom 
of Fig. 2(b)). Thus, we can describe a foliated triangula- 
tion by 'comparing' it to a reference lattice and specifying 
what forks are present and what forks are absent. This 
renders a description which is in spirit similar to that of 
a lattice gas model, where the description of a fluid, with 
molecules absent or present, is mapped to the description 
of a magnet with two-level spins on a fixed lattice [2]. 

Let us be more precise. Consider a 2D square lattice 
where a binary variable X nm is associated with every ver- 
tex (n, to) of the lattice, with 1 < n < N and 1 < to < M 
(see Fig. 2(a)). This forms a bit array A. To transform 
this bit array to a foliated triangulation, we put a fork 
on every site (n, to) it X nm = 1, and no fork if X nm = 0. 
This collection of forks defines a graph T = (V,E) in a 
natural way: each center of a fork defines a vertex v € V, 
and its legs become edges e e E of T. The space-like edge 
and the time-like edge pointing upwards connect to the 
first vertex to the left, and the time-like edge pointing 
downwards connects to the first vertex directly below or 
to the left (see Fig. 2). This can be formally described 
as follows. Let p be a rank 3 tensor with components 

Pn,m,m> ■= X nm > X [(to - m' (mod M)) + 1] . (5) 

Then the vertex at site (n,m) is adjacent to the vertices 
at sites (n+1, a), (n, 6), (n— 1, c) where a, b, c are defined 
by the inequalities 

< 

Pn+l,m—l,a < Pn+l,m—l,m' Vto' 6 [M], TO,' +^ a 
< Pn,m-l,b < Pn,m-l,m' Vto' € [M], m! ^ b 
< Pn-l,m,c < Pn-l,m,m' Vto' E [M], m' + 1 C, 

(6) 

where [M] = {1,2,..., M}. 

From triangulations to bit arrays. We now show how 
to map a foliated triangulation to a bit array. The recipe 
consists of six steps. First, construct the dual graph 
of T, denoted T* , whose basic building block is the d- 
fork (short for dual fork) . This consists of one time- like 
edge and two space-like edges "pointing to the right" (see 
Fig. 3(a,b)). Second, note that every d-fork crosses one 
space-like edge of T. We use this fact to attach the label 
Si € {1,2,..., Af} to a fork i if it crosses a space-like edge 
at time ts t (see Fig. 3(c), where the labels are also repre- 
sented as colors). Third, record the labels of the d-forks 
that appear in T* from left to right and write them down 
in an integer string S := (Si, S2, ■ ■ ■ , Sf) (cf. Diagram 
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10100001100100111110010010011010110000011010100110110111 

{0,11 

FIG. 2: From bit arrays to triangulations. The bit array (a) 
represents the foliated triangulation (b). Each bit stands for 
the presence / absence of a fork, which consists of one ver- 
tex, one space-like edge and two time-like edges. While the 
bits are fixed on a two-dimensional grid, their configuration 
determine the geometry, similar to a lattice gas model. (For 
illustration issues we have chosen non-periodic boundary con- 
ditions in both time and space.) To store the triangulation (b) 
on a double precision floating-point format (64-bit) machine 
a single data unit is required (c). 



(1)). This string representation has certain symmetries. 
Two integer strings correspond to the same dual triangu- 
lation if they are related by a sequence of the following 
operations: 

(i) commutations of contiguous integers if the inte- 
gers differ by at least two, i.e. (. . . , Si, ■ • •) ~ 
(...,S i+1 ,Si,...)ti \Si-S i+1 \>2; 

(ii) cyclic permutations, i.e. (S%, S2, ■ ■ ■ , Sp) ~ 
(Sp, Si, S2, ■ ■ ■ , Sp-x); 

(hi) inversion operation, i.e. (Si, S2, ■ ■ ■ , Sp) ~ 
(Sp, Sp-i, . . . , Si). 

Notice that (ii) and (iii) are only applicable for periodic 
boundary conditions on the spatial slices, while (i) is a 
degeneracy introduced by the integer string encoding. 

Fourth, define a string £ flat := (S N , S* N , . . . , S N ), 
consisting of complete integer-sequences, S* N := 
(1,2, ... , N). Apply operation (i) to arrange any S in 
successive integer-sequences, by allowing for incomplete 
sequences, e.g. (1, 3, 4, 6) (see Fig. 3(d)). Fifth, define the 
extended string S E by adding zeros to S wherever there 
is a missing element, e.g. (1,0,3,0,0,6) (see Fig. 3(e)). 
Generally, 



S := (^0,^,0^0,52,. 
r\ r 2 



(7) 



where A := (N, Si, S2, ■ ■ ■ , Sp, 1), and T is a vector with 
components Ti :— A i+ i — A^ + N — 1 (mod N), for i e 
{1,2,...,F+1}. 

Sixth, map S* E to A by arranging the entries of S in a 
two dimensional array of size N x M (going from bottom 



to top and left to right), and then replacing each positive 
entry by a 1. Formally, 



X l3 =e(Sf N _ l+ i) for 



1 < i < N 



(8) 



where Q(x) — if x — and Q(x) = 1 for x > (see 
Fig. 3(f)). The minimum lattice size necessary to encode 
all triangulations with at most to F forks with N spatial 
slice is given by N x M, where M = F — N + 1 (see 
Supplementary Information [17]). 

Note that by construction the integer encoding repre- 
sents a coordinate- free encoding of triangulations. Con- 
sequently, operations (i)-(iii) introduce an equivalence 
class on the space of all integer strings, and there exists 
a straightforward algorithm to single out its representa- 
tives. This can be illustrated with the help of a simple 
example, evaluating all unique histories for triangulations 
with 3 spatial slices from 3 to 5 forks. Starting from all 
possible strings, one can apply operations (i)-(iii) to sin- 
gle out one representative S of every equivalence class: 
(1,2,3), (1,2,3,3), (1,2,2,3), (1,2,3,2), (1,1,2,3), 
(1,2,3,3,3), (1,2,2,3,3), (1,2,3,2,3), (1,2,3,3,2), 
(1,2,2,2,3), (1,2,2,3,2), (1,2,3,2,3), (1,1,2,3,3), 
(1,1,2,2,3), (1,1,2,3,2), (1,2,1,2,3), and (1,1,1,2,3). 
See Supplementary Information [17] for a more thorough 
discussion. 

Applications of binary encoding. We now rephrase var- 
ious quantities of CDT in binary language. Consider 
first the Ricci scalar (Eq. (2)). The coordination number 
of a vertex at site (n, m) is 3 (because each fork has 3 
edges) plus additional A nm edges which depend on the 
surrounding bit array; explicitly: 



F(d-m) 

E 



[Ki+l,F(m+j-l) 

)] (9) 



where F(m!) := m! (mod M), and d is the index of the 
first non-zero entry to the left of m given by < p(n, m + 
1, d) < p(n, m + 1, rri) VW ^ d (see Eqs. (5), (6)). The 
Ricci scalar R nm is then given by 



Ft,, 



3- A r , 
3 + A r , 



(10) 



The action (Eq. (3)) is also easily computed in the 
binary encoding. To compute the Euler characteristic, 
note that every fork is associated to one vertex, 3 edges 
and 2 faces, hence it does not contribute to the value 
of x- Thus only forks relative to the boundary (that is, 
those which are placed at the boundary, either in the 
space or time dimension) contribute to x- For instance, 
for a torus (topology S 1 x S 1 ), x = 0. The volume is 
then given by twice the number of forks, N2 — 2F. This, 
together with C(T), which can be evaluated exactly for 
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FIG. 3: From triangulations to bit arrays. The triangulation T of panel (a) (which is the same as in Fig. 1(b)) is mapped 
to the bit array of panel (g). (a) We first construct the dual of T (a), and redraw it to give it the appearance of a brick wall 
(b), where every vertical (colored) line represents a dual fork (d-fork). The dual graph is mapped to a string S (c). Using 
operations (i)-(iii) we reorder S in successive incomplete integer-sequences (d), and extend the string by adding zeros where 
the integer-sequence is incomplete (e). Finally, we map this extended string to the bit array with Eq. (8) (f,g). (Notice that 
we consider here open boundary conditions on the spatial slices.) 



1+1 dimensional triangulations [15], suffices to evaluate 
the discretized path integral (Eq. (4)). 

Finally, the Pachner moves can also be expressed in the 
integer encoding (see Supplementary Information [17]). 
This allows us to show that the Pachner moves in 1+1 
dimensions are ergodic, i.e they generate all simplicial 
triangulations. 

Conclusion and Outlook. We have presented a binary 
encoding of discrete space-time in 1+1 dimensions. To 
this end, we have borrowed the discretization method of 



CDT and shown how to compute various quantities of in- 
terest of this theory in the binary description. Our results 
enable us to express a classical theory of discrete space- 
time as a lattice gas model. This approach has several 
potential applications. To begin with, it provides a nat- 
ural framework for quantization. Forks, the elementary 
geometrical unit, constitute independent classical degrees 
of freedom or 'normal modes' of discrete space-time. In 
the quantized theory, the presence or absence of a fork 
will be interpreted as elementary excitations of a quan- 
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turn lattice gas. These excitations will give rise to quan- 
tum fluctuations and, more generally, to quantum states 
representing superpositions of different space-times. De- 
spite using a time foliation, our approach to quantizing 
space-time is conceptually different from CDT quantum 
gravity, which is based on a path integral formulation. 

Discretization of space-time was introduced as a regu- 
larization tool to compute the path integral of general rel- 
ativity. If one considers discrete space-time as real, then 
forks could be seen as fundamental constituents of space- 
time, namely as basic geometric structures that connect 
events. A quantum formulation of this theory would nat- 
urally introduce quantum correlations in space-time. 

Finally, from an experimental perspective, by imple- 
menting quantum lattice gas models in the way we have 
introduced them, one could realize quantum simulations 
of fluctuating space-time. Such quantum simulators, ana- 
log or digital, could be used to simulate models of quan- 
tum gravity in modern quantum optics laboratories, e.g. 
using ultra-cold atoms in optical lattices or ion-traps [16]. 
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We present supporting material for our paper. The focus is on specific details of, and relationships between, the 
various encodings of triangulations: (I) we discuss the efficiency of the binary encoding, and derive the minimum 
lattice size necessary to encode all triangulations with at most to F forks distributed over N spatial slices; (II) we will 
translate the Pachner moves to operations on the integer string encoding. The Pachner moves are transformations in 
the set of triangulations which are ergodic, that is, any two triangulations are related by a finite sequence of Pachner 
moves. This connection can be used to verify the ergodicity of the Pachner moves in 1+1 dimensions; (III) we show 
how to utilize the integer string encoding and the corresponding symmetry operations introduced in the main paper 
to single out integer strings representing unique triangulations; and (IV) we compare the degeneracies arising from 
the binary and integer string encoding. 



Lattice size in binary encoding 



The minimum lattice size necessary to encode all triangulations T up to F forks distributed on N spatial slices is 
given by N x M , where M = F — N + 1. To show this, let us temporarily put aside the requirement that there is at 
least one fork on every spatial slice, equivalent to the requirement that S contains every element in {1,2,..., N} at 
least once. It is easy to see that S = {X, X, . . . , X}, where all elements X £ {1,2,..., N} are the same, maximizes 
r = (N-X,N-1,N-1,...,N-1,X-1), where T is a vector of length F + l, such that F E = F + (N - 1)F = FN 
and M — F. Consequently, by reintroducing the constraint to have at least one fork per spatial slice, the string S 
that will result in the largest M, has as many elements of the same kind, that are F — N + 1 d-forks on a single spatial 
slice X. It is best to study this case using the dual brick picture, where it is straightforward to see that d-forks, 
X ± n, for n > 1, can be moved freely in both directions by using operations (i)-(iii) (see Fig. 4). Therefore any of 
the N — 1 d-forks can be absorbed in a sequence that contains already one X, and thus the number of sequences 
necessary to encode all possible triangulations with F forks distributed on N spatial slices is F — N + 1. According 
to the mapping, see Eq. (7) in the main text, every sequence in 5 E corresponds to one row in the corresponding A. 



=F =E =E= ^ 

| — -] | — -3nni — ^nmr 



rz =E 3 



(a) (b) (c) (d) 

FIG. 4: Four different dual triangulating are represented (a-d). It is shown, that by successive application of operations (i)-(iii) 
of the main article acting on (a) one obtains (d) . The the d-forks colored in red are moved with respect to the previous dual 
triangulations. 



Integer equivalent of Pachner moves 



The Pachner moves are certain operations that when applied on a triangulation T generate other triangulations T". 
They are claimed to be ergodic, i.e. for any pair of triangulations, there always exists a sequence of Pachner moves 
that one can apply on one triangulation that yield the other. For foliated triangulations in 1+1 dimensions there are 
only two such moves (called Rule 1 and Rule 2, see Fig. 1 of the main text). Here we define the Pachner moves in the 
integer string description, that is, transformations on the integer strings which correspond to applying the Pachner 
move on the corresponding triangulation, as the following diagram illustrates. 



T 



Pachner 



V 



rp* 



(T'Y 



integer Pachner 1 ' 

5" 



(11) 
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(a) 





PACHNER MOVES 



S=(1 ,2,3,4,5,6,7,1 ,2,3,4,5,6,7,1 ,2,3,4,5,6,7,1 ,2,3,4,5,6,7,1 ,2,3,4,5,6,7,1 ,2,3,4,5,6,7,1 ,2,3,4,5,6,7) 



S'=(1 ,2,3,4,5,6,7,1 ,2,4,3,5,1 ,2,1 ,2,3,4,4,5,6,7,6,7,3,4,6,7,1 ,2,3,4,5,4,5,6,7,1 ,3,2,1 ,2,3,4,5,5,6,7,4,5,6,7) 



INTEGER ENCODING OF T 



INTEGER ANALOGUE OF PACHNER MOVES 



(b) 



,^S=(1234567 1234567 1234567 1234567 1234567 

' (1234567 12345 123465767 1234567 1234567 

(1234567 12345 123456767 1234567 1234567 

(1234567 12435 123456767 1234567 1234567 

(1234567 12435 12344567671234567 1234567 

(1234567 12435 12 132445676734567 1234567 

(1234567 12435 12 123445676734567 1234567 

(1234567 12435 12 123445676734567 1234567 

(1234567 12435 12 123445676734567 1234567 

(1234567 12435 12 123445676734567 1234567 

(1234567 12435 12 123445676734567 1234567 

(1234567 12435 12 123445676734567 12344567 

i (1234567 12435 12 12344567673467 123544567 

!— S'=(l 234567 12435 12 12344567673467 123454567 



1234567 1234567) 
1234567 1234567) 
1234567 1234567) 
1234567 1234567) 
1234567 1234567) 
1234567 1234567) 
1234567 1234567) 
1324MZ1g34567) 
132 flS5§fF74567) 
132 12345674567) 
132 1 234556 7 4 5 6 7) 
132 1 234556 7 4 5 6 7) 
132 1 234556 7 4 5 6 7) 
132 1 2 34556 7 4 5 6 7) 



INTEGER ENCODING OF T 



FIG. 5: We show explicitly how to apply successively Rl and R2 on T to obtain T 1 in both, the simplicial (a) and the 
integer string (b) encoding. In order to apply Rl it is sometimes necessary to use the pairwise commutation relations (PC), 
i.e. (. . . , Si, Si+i, ■ ■ ■ ) ~ (. . . , Si+i, Si, . . .) if \Si — Si+i\ > 2. Notice that the integer strings S and S' correspond to the 
triangulations T and T 1 of Fig. 1 in the main text, respectively. 



The action of Rule 1 (Rl) is easily understood in the dual triangulation, where it simply exchanges the order of 
two neighboring consecutive d-forks. In the integer description it corresponds to the following operation 

(. . . , Si, S i+ i, . . . ) > (. . . , S i+ i,Si, . . .). (12) 

where \Si — Sj+il = 1. Rule 2 (R2) has an equally simple equivalent, 

(. . . , (n - l),n, (n + 1) . . . ) ^ (. . . , (n - 1), n, n, (n + 1), . . . ), (13) 

duplicating a d-fork and placing it right next to it. It is thus clear that Rl and R2 connect different equivalence classes 
in the string encoding. This demonstrates the ergodicity of the Pachner moves in 1+1 dimensions: starting from the 
minimal string S = (1,2,3, . . . ,N — 2,N — 1, N), successive application of Rl and R2 generate all possible strings, and 
hence all possible triangulations. Fig. 5 shows the action of the "integer Pachner moves" on triangulations mentioned 
in the main text. 



Representatives of equivalent classes in all three encodings 

To evaluate the path integral of CDT in any dimensions analytically, one has to first single out triangulations which 
are representatives of their equivalence classes. We demonstrate how this procedure can be done in integer string 
encoding. First, note that to encode all triangulations of size V = 2F distributed over N spatial slices, one has to 
generate all strings containing at least N forks and at most F of them containing each element of {1 . . . N} at least 
once. We then apply symmetry operations (i)-(iii) as illustrated in the following example. 

Consider all triangulations containing up to 5 forks distributed over 3 spatial slices. We require at least one fork 
on every spatial slice; the smallest (in terms of 2-volume) triangulations contain 3 forks. The procedure to find the 
representatives for the equivalent classes is the following: 



1. generate all possible integer strings of lengths 3, 4 and 5, with elements in Si £ {1, 2, 3}; 



2. consider only strings that contain every element {1, 2, 3} at least once (so that only connected triangulations 
are taken into account); 



3. apply symmetry operations to single out unique triangulations / integer strings: if two integer strings are related 
by a sequence of the symmetry operations (i-ii) they are in the same equivalence class. One representative of 
each equivalence class give all unique histories / triangulations. It is best to demonstrate this at hand of an 
example. What are the representatives of triangulations containing 4 forks distributed over 3 spatial slices? 
From basic combinatorics we are left with nine valid integer string encodings: Si = (1, 2, 3, 3), S2 = (1, 3, 2, 3), 
S 3 = (1,3,3,2), S 4 = (1,2,2,3), S 5 = (1,2,3,2), S 6 = (1,3,2,2), S 7 = (1,1,2,3), S s = (1,1,3,2), and S 9 = 
(1, 2, 1, 3). Let us to begin with focus on the dual triangulations with one extra d-fork on the third row, i.e. S±, 
S2 and S3. Utilizing the symmetry operations we can show that S2 = (1, 3, 2, 3) = (3, 1, 2, 3) = (1, 2, 3, 3) = Si, 
and S 3 = (1,3,3,2) = (3,1,3,2) = (3,3,1,2) = (3,1,2,3) = (1,2,3,3) = Si, and consequently S 2 and S 3 
are redundant, and Si is a representative. Similarly we get Sg = (1,3,2,2) = (3,1,2,2) = (1,2,2,3) = S4, 
5 8 = (1, 1, 3, 2) = (1, 3, 1, 2) = (3, 1, 1, 2) = (1, 1, 2, 3) = S 7 , and S 9 = (1, 2, 3, 1) = (1, 1, 2, 3) = S 7 . Altogether 
we can for example single out the following representatives: Si, £4, S§, and S 7 . 

Next we present the binary and the integer description for this set of triangulations explicitly. For the binary 
description we need to consider arrays of size N x F — iV + 1 = 3 x 3 (see above). 



Representatives for 3 forks distributed over 3 spatial slices 
There is only one way to distribute 3 forks on 3 spatial slices, such that 

/ 1 

Si = (1,2, 3) Sf = (1,2,3, 0,0,0, 0,0,0) -> \(Sf ) = 10 0), see FIG. 6(1) (14) 

\ 1 



Representatives for 4 forks distributed over 3 spatial slices 
After applying the symmetry operations we are left with four representatives: 



S 2 = (1,2,3,3) Sf = (1,2,3, 0,0,3, 0,0,0) A(Sf ) = 10 0, see FIG. 6(2); 



5 3 = (1,2,2,3) Sf = (1,2,0, 0,2,3, 0,0,0) -» X(Sf ) = I 1 1 J, see FIG. 6(3); 

V 1 / 

Ci-V) 

5 4 = (1,2,3,2) -» Sf = (1,2,3, 0,2,0, 0,0,0) -» X(Sf ) = f 110, see FIG. 6(4) 



S 5 = (M,2,3) Sf = (1, 2, 3, 1, 0, 0, 0, 0, 0) A(Sf ) = I 1 I, see FIG. 6(5) 
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Representatives for 5 forks distributed over 3 spatial slices 
After applying the symmetry operations we are left with 12 representatives: 

1 1 1 

S* 6 = (1,2,3,3,3) -> Si = (1,2,3, 0,0,3, 0,0,3) -> X(Sf) = ( 1 

5 7 = (1, 2, 2, 3, 3) S E = (1, 2, 0, 0, 2, 3, 0, 0, 3) 

5 8 = (1, 2, 3, 2, 3) -> 5 8 E = (1, 2, 3, 0, 2, 3, X, X, X 

5 9 E = (1,2,3, 0,0,3, 0,2,0) -> A(S 9 E ) = 



1 

1 1 
A(S E ) = ( 1 1 



S 9 = (1,2,3,3,2) 
5i = (1,2,2,2,3 
5 n = (1,2,2,3,2 
S 12 = (1,2,3,2,3 
Sia = (1,1,2,3,3 
S 14 = (l,l,2,2,3 
Sib = (1,1,2,3,2 

516 = (1,2, 1,2, 3 

5 17 = (l,l,l,2,3 



dE 
"1 n 



(1,2,0, 0,2,0, 0,2,3) -> A(Sf 



A(S E 



S E =(1,2,0, 0,2,3, 0,2,0) 
-► £ E =(1,2,3, 0,2,3, 0,0,0) 

S E =(1,0,0, 1,2,3, 0,0,3) 
-> 5 E = (1,0,0, 1,2,0, 0,2,3) 
-)■ S E =(1,0,0, 1,2,3, 0,2,0) 
-> S E =(1,2,0, 1,2,3, 0,0,0) - MS 

5 E =(1,0,0, 1,0,0, 1,2,3) 



\(S$) = 



12.' 



Hsf 3 ) 



Hsf 5 ) 



A(5 E ) 



1 

1 1 

1 1 

1 

1 1 

1 1 

1 

1 

1 1 1 

1 

1 

1 1 1 

1 

1 1 

1 1 

1 

1 1 

1 

1 1 

1 

1 1 

1 1 

1 

1 1 

1 1 

1 

1 1 
1 1 

1 

1 

1 1 1 



see FIG. 6(6); 
see FIG. 6(7); 
see FIG. 6(8); 
see FIG. 6(9); 
see FIG. 6(10); 
see FIG. 6(11); 
see FIG. 6(12); 
see FIG. 6(13); 
see FIG. 6(14); 
see FIG. 6(15); 
see FIG. 6(16); 
see FIG. 6(17); 



(16) 



Degeneracies of the various encodings 

The binary and the integer encoding differ in the size of their configuration spaces. To show this, let us for a 
moment relax the constraint of only considering connected triangulations to get a rough estimate of the size of the 
configuration space. For F N we can simplify the minimum lattice size required to encode all triangulations to 
N x F. Thus the number of triangulations in the binary encoding is growing oc (2 N ) F , while the size of the integer 
encoding grows oc N F . Since both encodings are ergodic (in the sense of containing all triangulations with F forks 
distributed over iV spatial slices at least once) the binary encoding is less efficient than the integer encoding. The 
integer encoding has the additional advantage that the degeneracies of the configuration space are related to algebraic 
symmetry operations acting on the strings (see above). 




FIG. 6: The representatives of the equivalence classes for all triangulations up to 5 forks distributed on 3 spatial slices are 
shown. Notice that because of the periodic boundary conditions on the spatial slices, it is advantageous to display three periods 
(indicated by the dotted vertical lines). 



11 



This can be demonstrated with the help of the following example. Consider the integer string S — 
(1,2,3, 1,2,3, 1,2,3), representing a triangulation with coordination number 6 everywhere except at the bound- 
aries. Even for this trivial example the mapping to a binary encoding is not unique, if the lattice size is larger then 
3x3. In Fig. 7 we have depicted 6 seemingly different encodings of S. Fig. 7(a-d) are simply different binary encodings 
of S, while Fig. 7(e-f) are mappings of S to A after having applied various symmetry operations on S. 



(1,2,3, 1,2,3, 1,2,3) 



(1,2,3, 1,2,3, 1,2,3) 





I * 


6 \ 


4 






5 \ 


4 


(a) 









(1,2,3, 1,2,3, 1,2,3) 



(d) 






5 \ 


5 


3 




6 \ 


6 , 


4 




5 \ 


5 


4 


(b) 









(1 ,2,1 ,3,2,3,1 ,2,3) — > (1 ,2,3, 1 ,2,3, 1 ,2,3) 



— I -5^- * 


6 \ 


4 




5 \ 






4 


(e) 









(1,2,3, 1,2,3, 1,2,3) 




+ 
+ 
+ 



(1 ,2,3,1 ,2,1 ,3,2,3) — > (1 ,2,3, 1 ,2,3, 1 ,2,3) 





6\y \ 




5 




< + 








(f) 









FIG. 7: Six different binary encodings of S = (1,2,3, 1,2,3, 1,2,3) in a 3 x 4 lattice are represented. The numbers at the 
vertices are the coordination numbers. 



